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The discovery of microRNAs (miRNAs) has led to a paradigm shift in our basic 
understanding of gene regulation. Competing endogenous RNAs (ceRNAs) are the recent 
entrants adding to the complexities of miRNA mediated gene regulation. ceRNAs are 
RNAs that share miRNA recognition elements (MREs) thereby regulating each other. It 
is apparent that miRNAs act as rheostats that fine-tune gene expression and maintain 
the functional balance of various gene networks. Thus MREs in coding and non-coding 
transcripts have evolved to become the crosstalk hubs of gene interactions, affecting 
the expression levels and activities of different ceRNAs. Decoding the crosstalk between 
MREs mediated by ceRNAs is critical to delineate the intricacies in gene regulation, and 
we have just begun to unravel this complexity. 
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INTRODUCTION 

In recent years, advances in genomic sequencing have led to sev- 
eral discoveries that can be considered exceptions to the central 
dogma of molecular biology. Much of the DNA that does not 
encode proteins has been shown to code for various types of 
functional RNAs that have important regulatory roles. The com- 
plexities in gene regulation increased logarithmically following 
the identification of microRNAs (miRNAs). miRNAs are 18-22 
nucleotide long, evolutionarily conserved single-stranded RNA 
molecules that play a significant role in maintaining cellular iden- 
tity and homeostasis (Bartel, 2004). miRNA mediated regulatory 
networks were identified as drivers in diverse disease conditions 
including cancers, where they function either as oncogenes or 
as tumor suppressors (Hofacker, 2007; Kapsimali et al, 2007; 
Moffett and Novina, 2007; Croce, 2009). This was followed by the 
discovery of several non-coding RNAs with regulatory functions 
in disease conditions (Rinn et al, 2007; Clark and Mattick, 2011; 
Sun et al., 2013). Several studies have demonstrated the com- 
plex mechanisms by which miRNAs regulate gene expression by 
interacting with multiple networks (Hammond, 2007; Yang et al., 
2007; Sarver et al, 2010). This is because a single miRNA can reg- 
ulate tens or hundreds of mRNA targets and, conversely, several 
miRNAs can regulate a single mRNA (Eulalio et al, 2008; Bartel, 
2009; Friedman et al, 2009; Ebert and Sharp, 2012). Depending 
on the degree of homology of the seed sequence (6-8 nucleotide 
long miRNA region that determines target specificity), miRNAs 
can either completely degrade the target mRNAs (perfect base 



Abbreviations: ceRNA, competing endogenous RNAs; ciRs, circular RNAs; lncR- 
NAs, long non-coding RNAs; miRNA, microRNA; MREs, microRNA recognition 
elements; UTRs, Untranslated region. 



pairing) or act as road-blocks to translation (imperfect base pair- 
ing) (Yekta et al, 2004; Eulalio et al, 2008; Bartel, 2009). The 
multifaceted role of miRNAs can be appreciated by the recently 
described decoy activity of miRNAs, signaling a paradigm shift. 
Perrotti's group observed that in leukemic blasts, miR-328 is 
downregulated in a BCR/ABL-dependent manner (Eiring et al., 
2010). In these cells, there is increased expression of the RNA 
binding protein hnRNP E2, which interacts with mRNA for a 
transcription factor CEBPA, inhibiting its translation. miR-328 
functions as a decoy by competing with CEBPA mRNA for the 
hnRNP E2 binding site, thereby preventing the translational inhi- 
bition of CEBPA mRNA. This function is in addition to its 
canonical function of suppressing PIM1 protein kinase transla- 
tion, but nonetheless is critical to the regulatory control of the 
blast cell oncogenesis. 

Gene expression can be regulated in cis by enhancer sequences 
or in trans by genes that encode transcription factors (that act 
as activators or repressors) or RNA binding proteins. Recently, 
several studies have demonstrated that both coding and non- 
coding RNA molecules can regulate gene expression in trans by 
acting as sponges of miRNAs (Karreth et al, 2011; Salmena et al., 
201 1; Tay et al., 2011; Karreth and Pandolfi, 2013; Su et al., 2013). 
Together these RNA molecules are called "competing endogenous 
RNAs" (ceRNAs), which constitute a major proportion of gene 
regulators. The identification of this level of gene regulation 
can explain the correlation between genome size and increase 
in species complexities. Although the number of protein-coding 
genes is comparable between lower organisms and humans, there 
is a proportional increase between the number of non-coding 
transcripts and the complexity among species (Mattick, 201 la,b). 
These ceRNAs share sequences recognized by the miRNAs called 
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microRNA recognition elements (MREs). The ratio of noncoding 
to coding sequences in humans is about 47:1 (Frith et al, 2005) 
with up to 97% of the human genome constituting non protein- 
coding DNA (Tisseur et al., 2011). Further discovery of this 
additional layer of complexity in gene regulatory networks opens 
up new avenues for targeting therapy. The ceRNA concept, which 
describes the cross talk between RNA molecules mediated by 
miRNA recognition elements (MREs), is illustrated in Figure 1. 

miRNA SPONGES: A HISTORICAL PERSPECTIVE 

miRNAs regulate their target genes by binding to complementary 
sequences in the 3' UTR. In order to understand the function 
of these miRNAs in biological systems, a general approach is to 
design miRNA sponges or target mimics that harbor several bind- 
ing sites for a miRNA of interest. This concept of competitive 
target inhibition of miRNAs within a cell was first demonstrated 
by Phil Sharp's laboratory (Ebert et al, 2007). By this approach, 
one can decrease the number of miRNAs that are available to 
interact with an mRNA of interest, thus generating a loss-of- 
function phenotype in cell culture. At the same time, miRNA 
sponges were also observed to occur naturally. 

The first endogenous miRNA sponge was discovered in plants 
and was found to moderate a miRNA mediated response to envi- 
ronmental stress (Franco-Zorrilla et al., 2007). The identified 
sponge was a non-coding RNA, Induced by Phosphate Starvation 
1 (IPS1) with a binding site for miR-399, a phosphate starvation- 
induced miRNA. However, miR-399 binding does not induce 
degradation of the IPS1 transcript due to mismatched nucleotides 



in the binding site, but rather results in sequestration of miR-399 
from other target transcripts. Thus, IPS1 can effectively function 
as a sponge inhibiting the number of miR-399 molecules available 
for regulating its target PH02 mRNA. The Arabidopsis thaliana 
pho2 mutant is a phosphate over-accumulator. This mutant carries 
a mutation in the PH02 gene, encoding a ubiquitin-conjugating 
enzyme (UBC24), that leads to a reduction in full-length tran- 
scripts. (Aung et al., 2006; Lin et al, 2008). This mutation causes 
symptoms related to inorganic phosphate (Pi) toxicity, as a con- 
sequence of enhanced Pi uptake, facilitated translocation of Pi 
from roots to shoots and impaired Pi remobilization from old to 
your leaves (Bari et al., 2006; Chitwood and Timmermans, 2007). 
Initially during phosphate starvation, PH02 is downregulated 
(Bari et al, 2006; Chitwood and Timmermans, 2007). In response 
to phosphate starvation, IPS1 RNAs are induced. The latter can 
sequester miR-399 resulting in stabilization and increased accumu- 
lation of PH02 and, concomitantly, in reduced shoot phosphate 
content (Franco-Zorrilla et al., 2007). The authors coined the 
term "target mimicry" to define this mechanism of inhibition of 
miRNA activity. 

This idea of target mimicry demonstrates unanticipated com- 
plexity in the network of RNA regulatory interactions and was the 
first natural example of self-regulation by RNAs. Subsequently, 
target mimicry has been exploited to determine the function 
of plant miRNA families using a high-throughput approach to 
generate a large-scale collection of knockdowns for Arabidopsis 
thaliana miRNA families (Todesco et al., 2010). This idea was fol- 
lowed by the discovery of several mRNAs that can competitively 
inhibit regulatory small RNAs in prokaryotes (Figueroa-Bossi 
et al, 2009; Mandin and Gottesman, 2009; Overgaard et al, 2009). 
In marmoset T cells transformed with primate virus Herpesvirus 
saimiri (HSV), downregulation of host miRNAs miR-16, -27, and 
-142-3p by viral non-coding transcripts called HSURs (H. saimiri 
U-rich RNAs) was observed (Cazalla et al., 2010). Additionally, 
protein expression of the miR-27 target gene FOXOl was observed 
to be upregulated in the presence of HSUR1 , although neither the 
cellular compartment in which this interaction between HSUR 
and miRNA occurs nor the mechanism by which HSUR induces 
miRNA regulation are understood (Ebert and Sharp, 2010). 
Recently, sponge activity of viral RNA was confirmed by infection 
of mouse cells with murine cytomegalovirus resulting in rapid 
downregulation of anti-viral cellular miR-27 (Libri et al, 2012; 
Marcinowski et al, 2012). The authors in this study speculate 
that RNA mediated miRNA degradation could be a more general 
strategy employed by viruses to manipulate host cells. Another 
example for miRNA sponge activity in cancer progression is that 
of a non-coding RNA called "highly upregulated in liver cancer" 
(HULC), one of the most up regulated of all genes in hepatocel- 
lular carcinoma (Panzitt et al., 2007). Subsequently, Wang et al. 
documented that CREB (cAMP reponse element binding) protein 
is involved in upregulation of HULC and that HULC sequesters 
several miRNAs including miR-372, leading to derepression of its 
target protein, PRKACB, which in turn induces phosphorylation 
and activation of CREB (Wang et al, 2010). Thus HULC is an 
autoregulatory long non-coding RNA. 

In spite of these observations of the occurrence of competitive 
inhibition, it was not until 2009 that a wider biological function 
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FIGURE 1 | Schematic representation of ceRNAs mediated gene 
regulation. The regulation of miRNAs and mRNAs are mutual; as a 
consequence, levels of mRNA(s) containing relevant MREs could affect the 
levels and activity of another mRNA. The cross talk between mRNAs 
mediated by the MREs depends on the number and frequency of specific 
type of MREs. This cross talk is significantly influenced by the number and 
frequency of relevant MREs in the transcriptome of a given cell. The 
ceRNA-miRNA-mRNA cross talk maintains the overall activity and 
functional balance of gene networks in a cell. The ceRNA mediated 
regulation of mRNA-Y is illustrated here. mRNA-Y is a functional target of 
miR-B. However, the expression levels of mRNA-X that contains MREs for 
miRNA-B can act like a sponge to sequester the effect of miR-B on 
mRNA-Y (indicated with bold red line). Hence a decrease in expression 
levels of mRNA-X can also decrease the levels of mRNA-Y. 
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was proposed for this regulatory mechanism. Following the 
discovery of natural miRNA sponges it was speculated that these 
regulatory RNAs exist to maintain periodicity in gene expression, 
by allowing brief spurts of miRNA activity followed by a pause 
during which time target mRNA levels can recover (Chitwood 
and Timmermans, 2007). Further, Seitz observed three paradoxes 
with respect to miRNA mediated gene regulation: ( 1 ) miRNA reg- 
ulates several gene targets, yet the repression of only a few target 
genes significantly affects physiological function; (2) the extent 
of miRNA-mediated repression is remarkably lower as compared 
to endogenous genetic variations, which are well-tolerated; and 
(3) miRNA gene targets are significantly conserved among related 
species, but vary greatly between more distant species (Seitz, 
2009). Seitz proposed that several of the computationally iden- 
tified targets of miRNA can function as its competitive inhibitor, 
sequestering that miRNA, thereby preventing it from binding to 
its authentic mRNA target (Seitz, 2009). Seitz study postulates 
that miRNAs bind not only to their legitimate targets, whose 
activity is sensitive to small reduction in protein expression, but 
also to its pseudotargets that, on the contrary, are insensitive to 
such small reductions that can buffer and dilute miRNA activity. 
This formed the basis for the ceRNA hypothesis (Salmena et al., 
2011), which was substantiated by the identification and charac- 
terization of several novel categories of ceRNAs in mammalian 
tissues. In the following section, we describe the various types of 
ceRNAs and provide a review of the experimental and/or theoret- 
ical evidences for their regulatory function. These studies demon- 
strate the competition between RNA molecules containing same 
MREs, leading to the sharp, rapid, and dynamic upregulation of 
target mRNAs. 

PROTEIN CODING RNAs as ceRNAs 

A flurry of recent studies in mammalian cells and tissues demon- 
strated the competition between protein coding targets of miR- 
NAs (Cesana et al., 2011; Karreth et al., 2011; Sumazin et al., 
2011; Tay et al., 2011). These studies support a regulatory role 
for mRNAs independent of their protein coding function. Tay 
et al. tested the hypothesis that RNAs actively regulate each other 
through competitive sequestration of miRNA binding (Tay et al., 
2011). The authors identified and validated several endogenous 
protein-coding transcripts that regulate phosphatase and tensin 
homolog (PTEN), using a combined computational and experi- 
mental approach. PTEN is a tumor suppressor gene, mutations 
in which have been implicated in several cancers. Further sub- 
tle changes in PTEN dose can dictate critical outcomes with 
respect to tumor initiation and progression (Trotman et al., 2003; 
Alimonti et al, 2010; Berger et al, 2011). Thus identifying ceR- 
NAs for PTEN is of functional significance. This combinatorial 
analysis identified over 100 protein coding mRNAs that share 
MREs with PTEN. Furthermore, the authors observed the com- 
peting transcripts such as CNOT6L and VAPA to antagonize 
PI3K/AKT signaling, thereby contributing to growth- and tumor- 
suppressive properties. One additional PTEN ceRNA that has 
been functionally validated is Zinc finger E-box-binding home- 
obox 2 (ZEB2) (Karreth et al., 201 1). Using Sleeping Beauty inser- 
tional mutagenesis screen in a mouse model of melanoma, the 
authors identified 33 putative PTEN ceRNAs. Further validation 



was performed using RNAi-mediated gene silencing in human 
melanoma cells. In this manner, silencing of 6 ceRNAs, including 
ZEB2, showed significant attenuation in PTEN levels in a miRNA- 
dependent, protein coding-independent manner. In addition to 
the regulatory function of the ZEB2 transcript, its protein func- 
tion is well established as a transcription factor. The ZEB2 
transcription factor is an activator of epithelial-to-mesenchymal 
transition (EMT) and thus it is oncogenic and involved in the pro- 
gression of epithelial cancers (Vandewalle et al., 2005; Gregory 
et al, 2008; Vandewalle et al, 2009). On the other hand, in this 
study the authors show that reduction of ZEB2 mRNA expression 
activates the PI3K/AKT pathway and enhances cell transforma- 
tion. Reduction of ZEB2 mRNA expression commonly occurs in 
human melanomas and other cancers expressing low PTEN levels. 
Thus the mRNA and protein of ZEB2 have contradictory tumor 
suppressive and oncogenic functions, respectively. However, in 
light of a recent study showing that ZEB2 protein is oncosuppres- 
sive in melanoma (Caramel et al, 2013), the functional role of 
ZEB2 is context and tumor type dependent. 

In a third study, large scale analysis of gene expression data 
with matched miRNA expression profiles in human glioblastomas 
identified ~7000 genes that can act as miRNA sponges or target 
decoys (Sumazin et al., 201 1). Further biochemical analyses in cell 
lines demonstrated an intricate network regulating established 
driver genes including PTEN, PDGFRA, RBI, VEGFA, STAT3, 
and RUNX1, that allows cross talk between several canonical 
oncogenic pathways by titrating the activity of common miRNAs. 
Thus, although individual miRNAs contribute to small changes in 
gene expression, the combinatorial effect of multiple ceRNAs that 
interact via common MREs can be substantial and can explain 
the gene regulatory process in situations where there are no direct 
transcriptional or post-transcriptional interactions. While, MREs 
in a given transcript can regulate their own levels (cis regulation at 
transcript level) it can be noted that these MREs can also influence 
the levels and activity of other mRNAs (trans regulation) by com- 
peting for the miRNAs. mRNA mediated cis and trans regulation 
is schematically represented in Figure 2. 

PSEUD0GENES as ceRNAs 

The identification of a novel biological role for the expressed 
pseudogene of PTEN was the first experimental proof for the 
cross talk between coding and non-coding RNAs (Poliseno et al., 
2010). Poliseno et al. showed that in addition to protein-coding 




FIGURE 2 I c/s and trans regulation of mRNA. The MREs in the 3'UTR of 
mRNA can act as cis regulatory elements affecting its own stability and 
levels via miRNAs. On the other hand, MREs in the mRNAs functioning as 
ceRNAs can regulate the levels and activity of another mRNA by trans 
regulation. 



www.frontiersin.org 



January 2014 | Volume 5 | Article 8 | 3 



Kartha and Subramanian 



Competing endogenous RNAs (ceRNAs) 



function, RNAs possess regulatory function and demonstrated 
this using mRNAs produced by the PTEN tumor suppressor gene 
and its pseudogene PTENP1 (also known as PTENpgl, PTEN2, 
PTENi|/l). They showed that MREs are conserved between PTEN 
and PTENP1 and, as a consequence, overexpression of PTENP1 
3'UTR correspondingly increased cellular levels of PTEN by 
sequestering all miRNAs that inhibit PTEN. This in turn led to 
growth inhibition. The authors also observed a direct correla- 
tion between the expression of PTEN and PTENP1 in normal 
human tissues and prostate tumor samples suggesting these to 
be co-regulated. Further there was a direct relationship between 
PTENP1 copy number and PTEN expression in colon cancer 
tissues, indicating that PTENP1 transcript levels can regulate 
PTEN expression and thus act as a tumor suppressor (Poliseno 
et al., 2010). A similar observation was also reported in this 
study for KRAS gene: KRAS1P pseudogene pairs. The tran- 
script levels of KRAS and KRAS1P were positively correlated 
in prostate cancer tissues. Interestingly, amplification of the 
KRAS IP locus is reported for several human tumors, including 
neuroblastoma, retinoblastoma, and hepatocellular carcinoma 
(Plantaz et al, 1997; Zimonjic et al, 1999; Van Der Wal et al, 
2003). These observations indicate a proto-oncogenic role for 
KRAS IP. 

Supporting the regulatory function of pseudogenes, Poliseno 
et al. also reported conservation in the MREs of specific miR- 
NAs between a gene and its pseudogene, in the majority of 
the genes analyzed. For instance, OCT4 pseudogenes OCT4- 
pgl, 3, 4, and 5 maintained validated binding sites for miR- 
145, CDK4 pseudogene CDK4PS for miR-34 family, FOX03 
pseudogene FOX03B for miR-182, KRAS pseudogene KRAS1P 
maintained MREs for miR-143 and let-7 family, while E2F3 
pseudogene E2F3P1 maintained validated binding sites for miR- 
17 family but not for miR-34 family. DNMT3A pseudogene 
DNMT3AP1 did not maintain MREs for miR-29 family and 
miR-143. Thus it seems that the function of a pseudogene gen- 
erally mirrors its kindred gene function and contributes to a 
miRNA decoy mechanism (Poliseno et al., 2010). Pseudogenes 
can possess many of the same MREs that are harbored on the 
ancestral genes and thus act as perfect sponges (Salmena et al., 
2011). It is noteworthy that many ribosomal genes have several 
pseudogenes that are differentially regulated (Balasubramanian 
et al., 2009), indicating intricately dynamic mechanisms for RNA 
cross talk. 

Recently Johnsson et al. reported a hitherto unknown mech- 
anism of regulation of PTEN transcription and translation by 
pseudogenes (Johnsson et al., 2013). The authors identified and 
characterized two antisense RNA (asRNA) isoforms, a and p 
which are encoded by pseudogene, PTENP1, in addition to the 
sense transcript that was shown to act as a miRNA sponge. 
The asRNA was generally expressed at higher copy numbers as 
compared to their sense counterpart (Johnsson et al, 2013). 
The a isoform of asRNA, acting in trans was shown to nega- 
tively regulate PTEN transcription by recruiting repressive chro- 
matin remodelers, DNA methyltransferase 3A (DNMT3a) and 
Enhancer of Zeste (EZH2) to PTEN promoter region and cat- 
alyzes the formation of H3K27me3 chromatin mark (Johnsson 
et al., 2013). The isoform p\ on the other hand, regulated PTEN 



at the post-transcriptional level. The f5 asRNA isoform interacted 
with the PTENP1 sense transcript in the cytoplasm maintaining 
its stability, thereby affecting the sense miRNA sponge activ- 
ity. This study is an elegant example of the intricate regulatory 
network involving pseudogenes and coding RNAs. 

Another pseudogene acting as a miRNA sponge is OCT4-pg4, 
which is abnormally activated in hepatocellular carcinoma (Wang 
et al., 2013a). OCT4 is a POU transcription factor which is a 
pleiotropic regulator of gene expression in embryonic stem cells. 
The authors demonstrate that this pseudogene can protect OCT4 
transcripts from miR-145 inhibition, thereby upregulating its 
protein expression. Further a significant correlation was observed 
between high OCT4-pg4 levels and poor prognosis, thus indicat- 
ing this pseudogene to exert an oncogenic role in hepatocellular 
carcinoma. 

LONG N0NC0DING RNAs (IncRNAs) as ceRNAs 

Any expressed RNA that is more than 200 nt in length and 
lacks obvious protein coding capacity is considered a IncRNA 
(Ponting et al., 2009). The dysregulation of IncRNAs has increas- 
ingly been linked to many human diseases, including cancers 
(Shi et al, 2013). The important roles played by IncRNAs in 
brain development, neuron function and maintenance, and neu- 
rodegenerative diseases are also becoming increasingly evident 
(Wu et al., 2013). Specific IncRNAs function to regulate gene 
expression either through epigenetic mechanisms or by post- 
transcriptional events such as mRNA processing and degra- 
dation by interacting with splicing factors and 3' UTR ele- 
ments, respectively (Guttman et al, 2009; Khalil et al., 2009; 
Huarte et al., 2010). Examples of IncRNAs involved in epige- 
netic regulation include X-inactive specific transcript (XIST) and 
large intergenic noncoding RNAs (lincRNAs). Notably, recent 
global analysis of Argonaute (Ago) bound transcripts through 
high-throughput sequencing of RNAs isolated by crosslinking 
immunoprecipitation (HITS-CLIP) technique suggests miRNA 
mediated regulation of IncRNAs (Licatalosi et al, 2008; Chi et al., 
2009). 

The ability of IncRNAs to function as ceRNA was first demon- 
strated in muscle differentiation (Cesana et al., 201 1). The authors 
identified a muscle-specific IncRNA, linc-MDl, which governs 
the time of muscle differentiation by acting as a ceRNA in mouse 
and human myoblasts. This IncRNA sequestered miRNAs, miR- 
133 and miR-135 that specifically target protein coding genes 
MAML1 and MEF2C, needed to activate the differentiation pro- 
gram. Recently lincRNAs have also been shown to function as 
a miRNA sponge in human embryonic stem cell self-renewal 
(Cheng and Lin, 2013). Wang et al. showed that linc-RoR (reg- 
ulator of reprogramming) functions as a miRNA sponge to post- 
transcriptionally regulate mRNAs of the core transcription factors 
for maintaining embryonic stem cell pluripotency (Wang et al., 
2013b). Linc-RoR competes with mRNAs for OCT4, NANOG and 
SOX2 for miR-145 binding. Consistent with its sponge effect, linc- 
RoR copy number was observed to be much higher than that of 
miR-145 (>100 vs. 10-20 copies/cell) in self-renewing embry- 
onic stem cells compared to differentiating embryonic stem cells 
(20 vs. >500 copies/ cell). Thus the sponge effect of linc-RoR can 
disappear once the embryonic stem cells start to differentiate. 
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FIGURE 3 | Interactions between different types of ceRNAs, miRNA, 
and mRNA. Coding RNAs and non-coding RNAs including pseudogenes, 
IncRNAs, and circular RNAs can function as ceRNAs. 



CIRCULAR RNAs (circRNAs or ciRs) FUNCTION as ceRNAs 

Occasionally RNAs that are covalently linked at the ends to form 
circular molecules have been described in literature. circRNAs 
were discovered in plants and were shown to be derived from 
viroids (Sanger et al, 1976). In lower organisms, circRNAs stem 
from self-splicing introns of pre-ribosomal RNA (Grabowski 
et al, 1981). In animals, many circRNAs overlap with the cod- 
ing genes and arise from head-to-tail splicing of one or more 
exons. The best-known circRNA in animals is the antisense to the 
mRNA transcribed from SRY (sex-determining region Y) locus 
and considered circular testis-determining RNA (Capel et al., 
1993). Nevertheless, circRNAs were considered rare until recently 
when several groups demonstrated their widespread occurrence 
within transcriptomes (Salzman et al., 2012; Wu et al., 2012; Jeck 
et al., 2013). It was only very recently that a ceRNA function was 
attributed to this highly prevalent class of conserved regulatory 
RNA molecules (Hansen et al, 2013; Memczak et al., 2013). 

Memczak et al. adopted a novel computational methodology 
to screen RNA-seq libraries and detected thousands of human, 
mouse and C. elegans circRNA candidates (Memczak et al., 2013). 
They observed sequence conservation in these circRNAs and 
found several of these molecules to be specifically expressed in 
a cell type or developmental stage. Thus the authors speculated 
these previously disregarded molecules to perform a biological 
function, possibly in mediating gene regulation. Their hypothesis 
was proved by the characterization of a human circRNA, anti- 
sense to the cerebellar degeneration-related protein 1 transcript 
(CDRlas) also known as ciRS-7 (Hansen et al., 2013; Memczak 
et al., 2013). This brain-enriched circRNA harbors >70 MREs for 
the ancient miR-7 and forms complexes with Argonaute (AGO) 
protein in a miR-7 dependent manner. CDRlas (ciRS-7) is co- 
expressed with miR-7 in the brain and is co-localized in the 
P-bodies indicating that miR-7 compartmentalizes this circRNA 
to this organelle (Hansen et al., 2013). Further using functional 
approaches, it was demonstrated that the main role of CDRlas 
ciRS-7 is to bind miR-7 in neuronal tissues and ectopic expression 
of CDR1 is associated with morphological defects in the mid- 
brain, phenocopying miR-7 knock-down (Hansen et al., 2013; 
Memczak et al, 2013). 

It is noteworthy that non-coding RNAs, including pseudo- 
genes, IncRNAs, and circRNAs behave more effectively as ceRNAs 
than do protein-coding RNAs since non-coding RNAs are not 
involved in the active protein translation process (Gu et al., 
2009). While the miRNA sponges or miRNA-competing tran- 
scripts (protein-coding RNAs, pseudogenes, and IncRNAs) are 
expressed at low levels, contain a limited number of MREs, and 
are themselves subject to miRNA-mediated destabilization, cir- 
cRNAs have been shown to be long-lived in vivo compared to 
their linear counterparts and are completely resistant to miRNA- 
mediated target destabilization (Hansen et al, 2013; Memczak 
et al., 2013). A recent study also showed that circRNAs can be 
used as "next generation" exogenous sponges that outperform 
linear ones (Liu et al., 2013a). Circular miRNA sponges against 
miR-21 or miR-221 obtained using the self-splicing permuted 
intron-exon sequences derived from T4 bacteriophage gene td 
(Liu et al., 2013b) were introduced into malignant melanoma 
cell lines and resulted in excellent anti-cancer effects. In sum, 



the non-coding RNAs vast outnumber the protein coding genes 
and are extremely versatile in their mechanism of action (Djebali 
et al, 2012; Guttman and Rinn, 2012). A schematic representation 
of various types of ceRNAs and their unifying role in tweaking 
mRNA expression level and activity is provided in Figure 3. 

IDENTIFICATION OF ceRNAs 
COMPUTATIONAL METHODS 

The knowledge of precise numbers, types, and positions of MREs 
are very important for identification of ceRNAs. Although sev- 
eral target prediction algorithms are available, many are not 
successful in identifying all the important targets as the criteria 
for identifying targets are still evolving (Bartel, 2009). Tay et al. 
developed an approach termed mutually targeted MRE enrich- 
ment (MuTaME), involving integrated computational analysis 
and experimental validation (Tay et al, 2011). The authors used 
Rna22 miRNA target prediction algorithm (Miranda et al., 2006) 
to generate the MuTaME scores for the protein coding transcripts 
in humans. The MuTaME program evaluates a potential ceRNA 
based on (a) the number of miRNAs that are shared between 
transcripts of interest, (b) the number of predicted MREs in the 
transcript for a particular miRNA and the sequence span in which 
they are represented, (c) the distribution of the MREs of a partic- 
ular miRNA on the transcript and (d) the relation between the 
total numbers of MREs predicted in a transcript to the total num- 
ber of miRNAs that recognize these MREs. Using the MuTaME 
prediction algorithm Tay et al. (2011) demonstrated the phys- 
ical association of PTEN-targeting miRNAs with PTEN ceRNA 
CNOT6L. 

In order to assess the range and potential of miRNA-mediated 
regulatory interactions with a tumorigenic role, a novel multi- 
variate analysis method called Hermes was developed (Sumazin 
et al, 2011). Hermes uses a large collection of mRNA and 
miRNA expression profiles obtained from same tumor samples 
to identify modulators of miRNA activity. Further, in order to 
predict the ceRNA network, the HERMES algorithm also consid- 
ers non-coding RNAs. The program uses mutual and conditional 
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mutual information to quantify the dependency of one variable 
over the other variable, which in this case is affected by poten- 
tial ceRNAs. Using this program, the authors were able to predict 
the sponge interactions among ceRNAs in a glioblastoma dataset 
(Sumazin et al, 2011). 

Our group subsequently developed a database named "com- 
peting endogenous RNA database" (ceRDB) that can help identify 
ceRNAs for a given mRNA. The program examines the co- 
occurrence of MREs in the mRNAs (3' UTRs) on a genome-wide 
basis to predict ceRNAs for a specific mRNA targeted by miR- 
NAs (Sarver and Subramanian, 2012). The ceRNA finder database 
is available through www.oncomir.umn.edu/cefinder. The pro- 
gram uses over 50,000 conserved human miRNA-mRNA target 
interactions obtained from Targetscan predictions with a matrix 
containing 153 miRNA families and 9448 target mRNAs. For each 
mRNA, the interactions are scored based on the sum of the total 
number of miRNA binding sites for a given miRNA. The interac- 
tion scores based on the MRE frequency is sorted to predict the 
top 50 potential ceRNAs for a gene of interest. 

More recently a mathematical mass-action model to determine 
the optimal conditions for ceRNA activity in vivo was described by 
Pandolfi's group (Ala et al., 2013). The program mathematically 
models the complex network environment in which ceRNAs can 
function with optimal activity. Furthermore, the authors found 
that ceRNAs are responsive to modulations in the levels of various 
transcription factors. Thus, the deregulation of one ceRNA can 
have dramatic effect on the integrated ceRNA and transcriptional 
network. Besides the mass-action model there are also several 
other model systems reported (Bosia et al., 2013; De Giorgio et al., 
2013; Figliuzzi et al, 2013; Noorbakhsh et al, 2013). Linc2GO, 
a lincRNA (long intergenic non-coding RNA) function annota- 
tion database was recently generated to predict lincRNA functions 
based on ceRNA hypothesis (Liu et al., 2013a). 

The critical aspects of ceRNA identification are co-expression 
of the ceRNA, expression of miRNAs that commonly target the 
ceRNAs and the physical association (in part determined by the 
binding site availability due to RNA secondary structures) of miR- 
NAs with the ceRNAs. Although various miRNA target prediction 
programs use algorithms to include seed sequence matches, evo- 
lutionary conservation of binding sites, and flanking secondary 
structures (Bartel, 2009), they predominantly focus on the 3'UTR 
regions to predict the physical association between miRNAs and 
target genes. Based on the recent findings that pseudogenes, lncR- 
NAs and circular RNAs, can potentially function as ceRNAs, it 
becomes difficult to precisely determine the complete spectrum of 
all ceRNAs for a given gene of interest, unless non-coding RNAs 
are included in the target prediction algorithms. 

EXPERIMENTAL METHODS 

Several biochemical techniques have also enabled the exper- 
imental identification and functional validation of ceRNAs. 
For instance, Argonaute high-throughput sequencing of RNAs 
isolated by crosslinking immunoprecipitation (Ago HITS- 
CLIP), photoactivatable-ribonucleoside-enhanced crosslinking 
and immunoprecipitation (PAR-CLIP) and RNA immunopre- 
cipitation are some platforms that, by providing a footprint of 
all the miRNAs physically bound to a given RNA, can enable 



identification of MREs and will provide insights on ceRNA reg- 
ulatory networks (Thomas et al, 2010). 

ceRNAs: WHAT NEXT? 

It is becoming increasingly evident that regulation of gene expres- 
sion through competition for miRNA binding is a general phe- 
nomenon. The understanding of gene regulation at this com- 
plexity has provided a functional explanation on how certain 
previously known genes with no direct interactions are inter- 
linked, as exemplified in the studies involving PTEN. But still, 
there is more to be unraveled and learned. 

It is evident that the relative production and turnover rates of 
miRNAs, their target RNAs and ceRNAs can dictate the extent and 
duration of gene regulation (Arvey et al., 2010; Ala et al., 2013). 
There should be significant changes in ceRNA expression in order 
to relieve the miRNA repression of target mRNAs. Similarly, the 
magnitude of expression of the sequestered miRNAs should not 
be very low or abundant, as both conditions can obliterate the 
competition. Another factor that contributes to efficient compe- 
tition is the number of miRNAs that a ceRNA can sponge. This 
in turn depends on the subcellular distribution and interaction 
of ceRNAs with RNA-binding proteins and ribosomes. Both ceR- 
NAs and miRNAs should occur in the same tissue or cell type 
or cellular compartment at a specific time. Yet another factor that 
needs to be kept in mind is that not all MREs on ceRNAs are equal 
when it comes to binding a specific miRNA. This will depend on 
the nucleotide composition of MREs. All these details need to be 
clearly worked out to better appreciate the role of ceRNAs. 

Another potential area for future research is in understanding 
the role of ceRNAs in disease pathophysiology. Several stud- 
ies have analyzed gene or miRNA expression profiles in disease 
conditions. With the emergence of ceRNA mechanisms, what is 
warranted is a comprehensive analysis of both gene and miRNA 
expression in the same pathological sample. This will allow for 
identification of novel pathways that are deregulated in disease 
condition and decipher novel interactions between relevant sig- 
naling pathways. The bioinformatics tools that have been devel- 
oped in recent years can contribute largely to this amount of 
global analyses. 

Additionally, disease characterization in general involves 
extensive sequence analysis of coding regions to identify novel 
disease-causing mutations. These can be due to deletions, ampli- 
fications or chromosomal translocations as evident in cancer. 
With the emergence of ceRNAs, the search for perturbations 
should go beyond coding regions to include UTRs and introns 
as these can cause changes in MREs, which in turn can affect 
the complex regulatory circuits. For instance, gross genomic loss 
and amplification can have dramatic consequences for ceRNAs 
present in the region (Salmena et al., 2011). Similarly, chromo- 
somal translocations can lead to UTR swaps (Stephens et al., 
2009). Aberrant alternative splicing and alternative cleavage and 
polyadenylation events leading to shortening of 3'UTRs is another 
common mechanism in cancer (Mayr and Bartel, 2009). 

In addition to massive genomic rearrangements, point muta- 
tions can lead to genetic disorders as observed with inborn errors 
of metabolism. It is important to note that even though point 
mutations abolish protein function, the mutant transcripts retain 
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complete ceRNA function. Alternatively, point mutations in the 
MREs can abolish potential ceRNA functions, thereby perturbing 
the network. Sometimes these mutations are often overlooked, 
since they are not part of the functional protein coding sequence. 
Thus it is increasingly becoming important to have a whole 
genome approach when it comes to understanding disease. Often 
there is no correlation between the phenotype and genotype, 
suggesting genetic modifiers to contribute to the disease. With 
increasing evidence for miRNA mediated cross-talk, ceRNAs of 
a given disease-causing gene can potentially act as the modifier. 

MRE mediated cross-talk between transcripts is expected to 
change our approaches in understanding disease mechanisms. 
For example, in microarray data analysis, traditionally only signif- 
icant fold-changes are considered while interpreting gene expres- 
sion patterns. Knowing that the number of MREs and its dis- 
tribution can alter gene expression levels and its activity, even 
an insignificant fold-change in large number of genes is rele- 
vant and can considerably contribute a number of MREs that can 
affect the balance and dynamics of miRNA-mRNA interactions, 
which in turn affect cellular functions. Further, it is also possible 
that the transcripts with longer 3'UTR may function as "master 
MRE" containing genes and splice variants of these transcripts 
can also affect cell function. Gene regulation mediated by ceR- 
NAs can also impact targeted therapy and therapeutic response. 
Since transcripts can independently function as ceRNAs, target- 
ing their resulting proteins by small molecule inhibitors can in 
turn affect downstream the transcripts stability and levels that 
may contribute to undesirable side effects on treatment outcomes 
by the protein inhibitors. Modulation of levels of ceRNAs may 
be a strategic alternative to directly targeting a protein of inter- 
est. Further complicating this scenario is the recent finding that 
miRNAs can function as agents of intercellular communication 
(Arroyo et al., 2011; Lou et al., 2012; Steer and Subramanian, 
2012). It is noteworthy that exosomes, which recently gained a 
lot of attention with regard to its ability to mediate intercellu- 
lar communication (Mittelbrunn et al., 2011), may harbor both 
miRNAs and ceRNAs. This suggests that the miRNA pool in a 
cell is in constant flux, making it more challenging to understand 
gene regulation. Hence, it is critical to understand the spatial and 
temporal relationships between the relevant ceRNAs, mRNA, and 
miRNAs. 

CLOSING REMARKS 

ceRNAs mediated gene regulation is an emerging area of investi- 
gation that will significantly increase our understanding of disease 
pathobiology, particularly cancer. In general, multiple genes and 
signaling pathways are deregulated in cancer. Many of these 
genes have a redundant function that obscures identification of 
therapeutically relevant target genes. In addition, contradictory 
functions of mRNAs (as ceRNA) and proteins of a gene further 
challenge the identification of relevant target genes. The trans 
regulations of ceRNAs add yet another layer of complexity to 
gene regulation that needs to be carefully investigated for key 
oncogenes and tumor suppressors implicated in various cancer 
types. 

Overall, the normal functioning of cells depends on the main- 
tenance of balance between complex regulatory systems. This is 



achieved by dynamic and fine adjustments and through the ability 
to rapidly return to normal threshold levels when imbalances are 
detected. If the system fails to balance beyond threshold levels, 
the cells manifest a disease condition. The challenges we face is 
in understanding how this balance is maintained in physiolog- 
ical conditions and conversely, how it is altered in disease. The 
recent ENCODE project results suggest that about 85% of the 
human genome is transcribed, the majority of which are IncRNAs 
(Djebali et al, 2012; Hangauer et al, 2013). With the recent dis- 
covery of many novel small and long non-coding RNA types with 
distinctive functions, we are just beginning a new chapter in gene 
regulation and disease biology with more "unknown unknowns? 

ACKNOWLEDGMENTS 

We acknowledge the grant funding support from the Department 
of Defense (W81XWH10-1-0556) and the American Cancer 
Society (RSG-13-381-01) to Subbaya Subramanian. Due to space 
restrictions we could not cite many significant contributions 
made in this important and rapidly progressing field. 

REFERENCES 

Ala, U., Karreth, E A., Bosia, C, Pagnani, A., Taulli, R., Leopold, V., et al. 

(2013). Integrated transcriptional and competitive endogenous RNA networks 

are cross-regulated in permissive molecular environments. Proc. Natl. Acad. Sci. 

U.S.A. 110,7154-7159. doi: 10.1073/pnas.l222509110 
Alimonti, A., Carracedo, A., Clohessy, J. G., Trotman, L. C, Nardella, C, Egia, A., 

et al. (2010). Subtle variations in Pten dose determine cancer susceptibility. Nat. 

Genet. 42, 454-458. doi: 10.1038/ng.556 
Arroyo, J. D., Chevillet, J. R., Kroh, E. M., Ruf, I. K., Pritchard, C. C, Gibson, D. E, 

et al. (2011). Argonaute2 complexes carry a population of circulating microR- 

NAs independent of vesicles in human plasma. Proc. Natl. Acad. Sci. U.S.A. 108, 

5003-5008. doi: 10.1073/pnas.l019055108 
Arvey, A., Larsson, E., Sander, C, Leslie, C. S., and Marks, D. S. (2010). Target 

mRNA abundance dilutes microRNA and siRNA activity. Mol. Syst. Biol. 6, 363. 

doi: 10.1038/msb.2010.24 
Aung, K., Lin, S. I., Wu, C. C, Huang, Y. X, Su, C. L., and Chiou, T. J. 

(2006). pho2, a phosphate overaccumulator, is caused by a nonsense muta- 
tion in a microRNA399 target gene. Plant Physiol. 141, 1000-1011. doi: 

10.1104/pp.l06.078063 
Balasubramanian, S., Zheng, D., Liu, Y. J., Fang, G, Frankish, A., Carriero, N., et al. 

(2009). Comparative analysis of processed ribosomal protein pseudogenes in 

four mammalian genomes. Genome Biol. 10, R2. doi: 10.1186/gb-2009-10-l-r2 
Bari, R., Datt Pant, B., Stitt, M., and Scheible, W. R. (2006). PH02, microRNA399, 

and PHR1 define a phosphate-signaling pathway in plants. Plant Physiol. 141, 

988-999. doi: 10.1104/pp.l06.079707 
Bartel, D. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function. 

Cell 116, 281-297. doi: 10.1016/S0092-8674(04)00045-5 
Bartel, D. P. (2009). MicroRNAs: target recognition and regulatory functions. Cell 

136,215-233. doi: 10.1016/j.cell.2009.01.002 
Berger, A. EL, Knudson, A. G, and Pandolfi, P. P. (2011). A continuum model for 

tumour suppression. Nature 476, 163-169. doi: 10.1038/naturel0275 
Bosia, C, Pagnani, A., and Zecchina, R. (2013). Modelling Competing Endogenous 

RNA Networks. PLoS ONE 8:e66609. doi: 10.1371/journal.pone.0066609 
Capel, B., Swain, A., Nicolis, S., Hacker, A., Walter, M., Koopman, P., et al. (1993). 

Circular transcripts of the testis-determining gene Sry in adult mouse testis. Cell 

73, 1019-1030. doi: 10.1016/0092-8674(93)90279-Y 
Caramel, J., Papadogeorgakis, E., Hill, L., Browne, G. J., Richard, G, Wierinckx, 

A., et al. (2013). A switch in the expression of embryonic EMT-inducers 

drives the development of malignant melanoma. Cancer Cell 24, 466-480. doi: 

10.1016/j.ccr.2013.08.018 
Cazalla, D., Yario, T., and Steitz, J. A. (2010). Down-regulation of a host 

microRNA by a Herpesvirus saimiri noncoding RNA. Science 328, 1563-1566. 

doi: 10.1126/science.ll87197 
Cesana, M., Cacchiarelli, D., Legnini, I., Santini, T., Sthandier, O., Chinappi, 

M., et al. (2011). A long noncoding RNA controls muscle differentiation 



www.frontiersin.org 



January 2014 | Volume 5 | Article 8 | 7 



Kartha and Subramanian 



Competing endogenous RNAs (ceRNAs) 



by functioning as a competing endogenous RNA. Cell 147, 358-369. doi: 

10.1016/j.cell.201 1.09.028 
Cheng, E. C, and Lin, H. (2013). Repressing the repressor: a lincRNA as a 

MicroRNA sponge in embryonic stem cell self-renewal. Dev. Cell 25, 1-2. doi: 

10.1016/j.devcel.2013.03.020 
Chi, S. W., Zang, J. B., Mele, A., and Darnell, R. B. (2009). Argonaute HITS- 
CLIP decodes microRNA-mRNA interaction maps. Nature 460, 479-486. doi: 

10.1038/nature08170 
Chitwood, D. H., and Timmermans, M. C. (2007). Target mimics modulate 

miRNAs. Nat. Genet. 39, 935-936. doi: 10.1038/ng0807-935 
Clark, M. B., and Mattick, J. S. (20 1 1 ) . Long noncoding RNAs in cell biology. Semin. 

CellDev. Biol. 22, 366-376. doi: 10.1016/j.semcdb.2011.01.001 
Croce, C. M. (2009). Causes and consequences of microRNA dysregulation in 

cancer. Nat. Rev. Genet. 10, 704-714. doi: 10.1038/nrg2634 
De Giorgio, A., Krell, J., Harding, V., Stebbing, J., and Castellano, L. (2013). 

Emerging roles of competing endogenous RNAs in cancer: insights from the 

regulation of PTEN. Mol. Cell. Biol. 33, 3976-3982. doi: 10.1 128/MCB.00683-13 
Djebali, S., Davis, C. A., Merkel, A., Dobin, A., Lassmann, T, Mortazavi, A., et al. 

(2012). Landscape of transcription in human cells. Nature 489, 101-108. doi: 

10.1038/naturell233 
Ebert, M. S., Neilson, J. R., and Sharp, P. A. (2007). MicroRNA sponges: competitive 

inhibitors of small RNAs in mammalian cells. Nat. Methods 4, 721-726. doi: 

10.1038/nmethl079 

Ebert, M. S., and Sharp, P. A. (2010). Emerging roles for natural microRNA 
sponges. Curr. Biol. 20, R858-R861. doi: 10.1016/j.cub.2010.08.052 

Ebert, M. S., and Sharp, P. A. (2012). Roles for microRNAs in conferring robustness 
to biological processes. Cell 149, 515-524. doi: 10.1016/j.cell.2012.04.005 

Eiring, A. JVL, Harb, J. G., Neviani, P., Garton, C, Oaks, J. J., Spizzo, R-, et al. 
(2010). miR-328 functions as an RNA decoy to modulate hnRNP E2 reg- 
ulation of mRNA translation in leukemic blasts. Cell 140, 652-665. doi: 
10.1016/j.cell.2010.01.007 

Eulalio, A., Huntzinger, E., and Izaurralde, E. (2008). Getting to the root of miRNA- 
mediated gene silencing. Cell 132, 9-14. doi: 10.1016/j.cell.2007.12.024 

Figliuzzi, M., Marinari, E., and De Martino, A. (2013). MicroRNAs as a selective 
channel of communication between competing RNAs: a steady-state theory. 
Biophys.J. 104, 1203-1213. doi: 10.1016/j.bpj.2013.01.012 

Figueroa-Bossi, N., Valentini, M., Malleret, L., Fiorini, E, and Bossi, L. (2009). 
Caught at its own game: regulatory small RNA inactivated by an inducible tran- 
script mimicking its target. GenesDev. 23, 2004-2015. doi: 10.1101/gad.541609 

Franco-Zorrilla, J. M., Valli, A., Todesco, M., Mateos, I., Puga, M. I., Rubio-Somoza, 
I., et al. (2007). Target mimicry provides a new mechanism for regulation of 
microRNA activity. Nat. Genet. 39, 1033-1037. doi: 10.1038/ng2079 

Friedman, R. C, Farh, K. K., Burge, C. B., and Bartel, D. P. (2009). Most mam- 
malian mRNAs are conserved targets of microRNAs. Genome Res. 19, 92-105. 
doi: 10.1101/gr.082701.108 

Frith, M. C, Pheasant, M., and Mattick, J. S. (2005). The amazing complex- 
ity of the human transcriptome. Eur. J. Hum. Genet. 13, 894-897. doi: 
10.1038/sj.ejhg.5201459 

Grabowski, P. J., Zaug, A. J., and Cech, T. R. (1981). The intervening sequence of 
the ribosomal RNA precursor is converted to a circular RNA in isolated nuclei 
of Tetrahymena. Cell 23, 467-476. doi: 10.1016/0092-8674(81)90142-2 

Gregory, P. A., Bert, A. G, Paterson, E. L., Barry, S. C, Tsykin, A., Farshid, G-, 
et al. (2008). The miR-200 family and miR-205 regulate epithelial to mesenchy- 
mal transition by targeting ZEB1 and SIP1. Nat. Cell Biol. 10, 593-601. doi: 
10.1038/ncbl722 

Gu, S., Jin, L., Zhang, F., Sarnow, P., and Kay, M. A. (2009). Biological basis for 
restriction of microRNA targets to the 3' untranslated region in mammalian 
mRNAs. Nat. Struct. Mol. Biol. 16, 144-150. doi: 10.1038/nsmb.l552 

Guttman, M., Amit, I., Garber, M., French, C, Lin, M. E, Feldser, D., et al. (2009). 
Chromatin signature reveals over a thousand highly conserved large non-coding 
RNAs in mammals. Nature 458, 223-227. doi: 10.1038/nature07672 

Guttman, M., and Rinn, J. L. (2012). Modular regulatory principles of large non- 
coding RNAs. Nature 482, 339-346. doi: 10.1038/naturel0887 

Hammond, S. M. (2007). MicroRNAs as tumor suppressors. Nat. Genet. 39, 
582-583. doi: 10.1038/ng0507-582 

Hangauer, M. J., Vaughn, I. W., and McManus, M. T. (2013). Pervasive tran- 
scription of the human genome produces thousands of previously unidentified 
long intergenic noncoding RNAs. PLoS Genet. 9:el003569. doi: 10.1371/jour- 
nal.pgen. 1003569 



Hansen, T. B., Jensen, T. I., Clausen, B. H, Bramsen, J. B., Finsen, B., Damgaard, C. 

K., et al. (2013). Natural RNA circles function as efficient microRNA sponges. 

Nature 495, 384-388. doi: 10.1038/naturell993 
Hofacker, I. L. (2007). How microRNAs choose their targets. Nat. Genet. 39, 

1191-1192. doi: 10.1038/ngl007-1191 
Huarte, M., Guttman, M., Feldser, D., Garber, M., Koziol, M. J., Kenzelmann- 

Broz, D., et al. (2010). A large intergenic noncoding RNA induced by p53 

mediates global gene repression in the p53 response. Cell 142, 409-419. doi: 

10.1016/j.cell.2010.06.040 
Jeck, W. R., Sorrentino, J. A., Wang, K., Slevin, M. K., Burd, C. E., Liu, J., et al. 

(2013). Circular RNAs are abundant, conserved, and associated with ALU 

repeats. RNA 19, 141-157. doi: 10.1261/rna.035667.112 
Johnsson, P., Ackley, A., Vidarsdottir, L., Lui, W. O., Corcoran, M., Grander, D., 

et al. (2013). A pseudogene long-noncoding-RNA network regulates PTEN 

transcription and translation in human cells. Nat. Struct. Mol. Biol. 20, 440-446. 

doi: 10.1038/nsmb.2516 
Kapsimali, M., Kloosterman, W. P., De Bruijn, E., Rosa, E, Plasterk, R. H., and 

Wilson, S. W. (2007). MicroRNAs show a wide diversity of expression profiles 

in the developing and mature central nervous system. Genome Biol. 8, R173. doi: 

10.1186/gb-2007-8-8-rl73 
Karreth, F. A., and Pandolfi, P. P. (2013). ceRNA cross-talk in cancer: when ce-bling 

rivalries go awry. Cancer Discov. 3, 1113-1121. doi: 10.1158/2159-8290.CD-13- 

0202 

Karreth, F. A., Tay, Y., Perna, D., Ala, U., Tan, S. M., Rust, A. G, et al. (2011). 
In Vivo identification of tumor- suppressive PTEN ceRNAs in an onco- 
genic BRAF-induced mouse model of melanoma. Cell 147, 382-395. doi: 
10.1016/j.cell.201 1.09.032 

Khalil, A. M., Guttman, M., Huarte, M., Garber, M., Raj, A., Rivea Morales, D., 
et al. (2009). Many human large intergenic noncoding RNAs associate with 
chromatin-modifying complexes and affect gene expression. Proc. Natl. Acad. 
Sci. U.S.A. 106, 11667-11672. doi: 10.1073/pnas.0904715106 

Libri, V., Helwak, A., Miesen, P., Santhakumar, D., Borger, J. G-, Kudla, G, 
et al. (2012). Murine cytomegalovirus encodes a miR-27 inhibitor disguised 
as a target. Proc. Natl. Acad. Sci. U.S.A. 109, 279-284. doi: 10.1073/pnas.lll4 
204109 

Licatalosi, D. D, Mele, A., Fak, J. J., Ule, J., Kayikci, M., Chi, S. W., et al. (2008). 

HITS-CLIP yields genome-wide insights into brain alternative RNA processing. 

Nature 456, 464-469. doi: 10.1038/nature07488 
Lin, S. I., Chiang, S. E, Lin, W. Y., Chen, J. W., Tseng, C. Y, Wu, P. C, et al. (2008). 

Regulatory network of microRNA399 and PH02 by systemic signaling. Plant 

Physiol. 147, 732-746. doi: 10.1104/pp.l08.116269 
Liu, K., Yan, Z., Li, Y, and Sun, Z. (2013a). Linc2GO: a human LincRNA function 

annotation resource based on ceRNA hypothesis. Bioinformatics 29, 2221-2222. 

doi: 10.1093/bioinformatics/btt361 
Liu, Y, Cui, H., Wang, W., Li, L., Wang, Z., Yang, S., et al. (2013b). Construction 

of circular miRNA sponges targeting miR-21 or miR-221 and demonstration of 

their excellent anticancer effects on malignant melanoma cells. Int. }. Biochem. 

Cell Biol. 45, 2643-2650. doi: 10.1016/j.biocel.2013.09.003 
Lou, E., Fujisawa, S., Barlas, A., Romin, Y, Manova-Todorova, K., Moore, M. A., 

et al. (2012). Tunneling Nanotubes: a new paradigm for studying intercellular 

communication and therapeutics in cancer. Commun. Integr. Biol. 5, 399-403. 

doi: 10.4161/cib.20569 
Mandin, P., and Gottesman, S. (2009). Regulating the regulator: an RNA decoy acts 

as an OFF switch for the regulation of an sRNA. Genes Dev. 23, 1981-1985. doi: 

10.1 10 1/gad. 1846609 
Marcinowski, L., Tanguy, M., Krmpotic, A., Radle, B., Lisnic, V. J., Tuddenham, 

L., et al. (2012). Degradation of cellular mir-27 by a novel, highly abundant 

viral transcript is important for efficient virus replication in vivo. PLoS Pathog. 

8:el002510. doi: 10.1371/journal.ppat.l002510 
Mattick, J. S. (2011a). Genome-sequencing anniversary. The genomic foundation 

is shifting. Science 331, 874. doi: 10.1126/science.l203703 
Mattick, J. S. (2011b). Long noncoding RNAs in cell and developmental biology. 

Semin. CellDev. Biol. 22, 327. doi: 10.1016/j.semcdb.2011.05.002 
Mayr, C, and Bartel, D. P. (2009). Widespread shortening of 3'UTRs by alterna- 
tive cleavage and polyadenylation activates oncogenes in cancer cells. Cell 138, 

673-684. doi: 10.1016/j.cell.2009.06.016 
Memczak, S., Jens, M., Elefsinioti, A., Torti, E, Krueger, J., Rybak, A., et al. (2013). 

Circular RNAs are a large class of animal RNAs with regulatory potency. Nature 

495, 333-338. doi: 10.1038/naturell928 



Frontiers in Genetics | Non-Coding RNA 



January 2014 | Volume 5 | Article 8 | 8 



Kartha and Subramanian 



Competing endogenous RNAs (ceRNAs) 



Miranda, K. C, Huynh, T., Tay, Y., Ang, Y. S., Tam, W. L., Thomson, A. M., 
et al. (2006). A pattern-based method for the identification of MicroRNA bind- 
ing sites and their corresponding heteroduplexes. Cell 126, 1203-1217. doi: 
10.1016/j.cell.2006.07.031 

Mittelbrunn, M., Gutierrez- Vazquez, C, Villarroya-Beltri, C, Gonzalez, S., 
Sanchez-Cabo, E, Gonzalez, M. A., et al. (2011). Unidirectional transfer of 
microRNA-loaded exosomes from T cells to antigen-presenting cells. Nat. 
Commun. 2, 282. doi: 10.1038/ncommsl285 

Moffett, H. E, and Novina, C. D. (2007). A small microRNA makes a Bic difference. 
Genome Biol. 8, 221. doi: 10.1186/gb-2007-8-7-221 

Noorbakhsh, J., Lang, A. H., and Mehta, P. (2013). Intrinsic noise of microRNA- 
regulated genes and the ceRNA hypothesis. PLoS ONE 8:e72676. doi: 
10.1371/journal.pone.0072676 

Overgaard, M., Johansen, J., Moller-Jensen, J., and Valentin-Hansen, R (2009). 
Switching off small RNA regulation with trap-mRNA. Mol. Microbiol. 73, 
790-800. doi: 10.1111/j.l365-2958.2009.06807.x 

Panzitt, K., Tschernatsch, M. M., Guelly, C., Moustafa, T., Stradner, M., Strohmaier, 
H. M., et al. (2007). Characterization of HULC, a novel gene with striking up- 
regulation in hepatocellular carcinoma, as noncoding RNA. Gastroenterology 
132, 330-342. doi: 10.1053/j.gastro.2006.08.026 

Plantaz, D., Mohapatra, G., Matthay, K. K., Pellarin, M., Seeger, R. C, and 
Feuerstein, B. G. (1997). Gain of chromosome 17 is the most frequent abnor- 
mality detected in neuroblastoma by comparative genomic hybridization. Am. 
J. Pathol 150, 81-89. 

Poliseno, L., Salmena, L., Zhang, J., Carver, B., Haveman, W. J., and Pandolfi, 
P. P. (2010). A coding-independent function of gene and pseudogene 
mRNAs regulates tumour biology. Nature 465, 1033-1038. doi: 10.1038/nature 
09144 

Ponting, C. P., Oliver, P. L., and Reik, W. (2009). Evolution and functions of long 
noncoding RNAs. Cell 136, 629-641. doi: 10.1016/j.cell.2009.02.006 

Rinn, J. L., Kertesz, M., Wang, ). K., Squazzo, S. L., Xu, X., Brugmann, S. 
A., et al. (2007). Functional demarcation of active and silent chromatin 
domains in human HOX loci by noncoding RNAs. Cell 129, 1311-1323. doi: 
10.1016/j.cell.2007.05.022 

Salmena, L., Poliseno, L., Tay, Y, Kats, L., and Pandolfi, P. P. (2011). A ceRNA 
hypothesis: the Rosetta Stone of a hidden RNA language? Cell 146, 353-358. 
doi: 10.1016/j.cell.2011.07.014 

Salzman, J., Gawad, C, Wang, P. L., Lacayo, N., and Brown, P. O. (2012). Circular 
RNAs are the predominant transcript isoform from hundreds of human genes 
in diverse cell types. PLoS ONE 7:e30733. doi: 10.1371/journal.pone.0030733 

Sanger, H. L., Klotz, G., Riesner, D., Gross, H. J., and Kleinschmidt, A. K. (1976). 
Viroids are single-stranded covalently closed circular RNA molecules exist- 
ing as highly base-paired rod-like structures. Proc. Natl Acad. Sci. U.S.A. 73, 
3852-3856. doi: 10.1073/pnas.73.11.3852 

Sarver, A., Li, L., and Subramanian, S. (2010). MicroRNA miR-183 functions as an 
oncogene by targeting the transcription factor EGR1 and promoting tumor cell 
migration. Can. Res. 70, 9570-9580. doi: 10.1158/0008-5472.CAN-10-2074 

Sarver, A. L., and Subramanian, S. (2012). Competing endogenous RNA database. 
Bioinformation 8, 731-733. doi: 10.6026/97320630008731 

Seitz, H. (2009). Redefining microRNA targets. Curr. Biol 19, 870-873. doi: 
10.1016/j.cub.2009.03.059 

Shi, X., Sun, M., Liu, LL, Yao, Y, and Song, Y. (2013). Long non-coding RNAs: a 
new frontier in the study of human diseases. Cancer Lett. 339, 159-166. doi: 
10.1016/j.canlet.2013.06.013 

Steer, C. J., and Subramanian, S. (2012). Circulating microRNAs as biomarkers: a 
new frontier in diagnostics. Liver Transpl. 18, 265-269. doi: 10.1002/lt.23377 

Stephens, P. J., McBride, D. J., Lin, M. L., Varela, I., Pleasance, E. D., Simpson, J. T, 
et al. (2009). Complex landscapes of somatic rearrangement in human breast 
cancer genomes. Nature 462, 1005-1010. doi: 10.1038/nature08645 

Su, X., Xing, J., Wang, Z., Chen, L., Cui, M., and Jiang, B. (2013). microRNAs 
and ceRNAs: RNA networks in pathogenesis of cancer. Chin. J. Cancer Res. 25, 
235-239. doi: 10.3978/j.issn.l000-9604 

Sumazin, P., Yang, X., Chiu, H. S., Chung, W. J., Iyer, A., Llobet-Navas, D., et al. 
(2011). An extensive microRNA-mediated network of RNA- RNA interactions 
regulates established oncogenic pathways in glioblastoma. Cell 147, 370-381. 
doi: 10.1016/j.cell.2011.09.041 

Sun, L., Goff, L. A., Trapnell, C, Alexander, R., Lo, K. A., Hacisuleyman, E., et al. 
(2013). Long noncoding RNAs regulate adipogenesis. Proc. Natl. Acad. Sci. 
U.S.A. 110, 3387-3392. doi: 10.1073/pnas.l222643110 



Tay, Y., Kats, L., Salmena, L., Weiss, D., Tan, S. M., Ala, U, et al. (2011). 
Coding-independent regulation of the tumor suppressor PTEN by competing 
endogenous mRNAs. Cell 147, 344-357. doi: 10.1016/j.cell.201 1.09.029 

Thomas, M., Lieberman, J., and Lai, A. (2010). Desperately seeking microRNA 
targets. Nat. Struct. Mol. Biol. 17, 1169-1174. doi: 10.1038/nsmb.l921 

Tisseur, M., Kwapisz, M., and Morillon, A. (201 1). Pervasive transcription - Lessons 
from yeast. Biochimie 93, 1889-1896. doi: 10.1016/j.biochi.201 1.07.001 

Todesco, M., Rubio-Somoza, I., Paz-Ares, )., and Weigel, D. (2010). A collection of 
target mimics for comprehensive analysis of microRNA function in Arahidopsis 
thaliana. PLoS Genet. 6:el001031. doi: 10.1371/journal.pgen.l001031 

Trotman, L. C, Niki, M., Dotan, Z. A., Koutcher, J. A., Di Cristofano, A., Xiao, A., 
et al. (2003). Pten dose dictates cancer progression in the prostate. PLoS Biol. 
1:E59. doi: 10.1371/journal.pbio.0000059 

Van Der Wal, J. E., Hermsen, M. A., Gille, H. J., Schouten-Van Meeteren, N. Y, 
Moll, A. C, Imhof, S. M., et al. (2003). Comparative genomic hybridisation 
divides retinoblastomas into a high and a low level chromosomal instability 
group./. Clin. Pathol. 56, 26-30. doi: 10.1136/jcp.56.1.26 

Vandewalle, C, Comijn, J., De Craene, B., Vermassen, P., Bruyneel, E., Andersen, 
LL, et al. (2005). SIP1/ZEB2 induces EMT by repressing genes of dif- 
ferent epithelial cell-cell junctions. Nucleic Acids Res. 33, 6566-6578. doi: 
10.1093/nar/gki965 

Vandewalle, C, Van Roy, E, and Berx, G. (2009). The role of the ZEB family of tran- 
scription factors in development and disease. Cell. Mol. Life Sci. 66, 773-787. 
doi: 10.1007/s00018-008-8465-8 

Wang, J., Liu, X., Wu, LL, Ni, P., Gu, Z., Qiao, Y, et al. (2010). CREB up-regulates 
long non-coding RNA, HULC expression through interaction with microRNA- 
372 in liver cancer. Nucleic Acids Res. 38, 5366-5383. doi: 10.1093/nar/gkq285 

Wang, L., Guo, Z. Y, Zhang, R., Xin, B., Chen, R., Zhao, J., et al. (2013a). 
Pseudogene OCT4-pg4 functions as a natural micro RNA sponge to regu- 
late OCT4 expression by competing for miR-145 in hepatocellular carcinoma. 
Carcinogenesis 34, 1773-1781. doi: 10.1093/carcin/bgtl39 

Wang, Y, Xu, Z., Jiang, J., Xu, C., Kang, J., Xiao, L., et al. (2013b). 
Endogenous miRNA sponge lincRNA-RoR regulates Oct4, Nanog, and Sox2 
in human embryonic stem cell self-renewal. Dev. Cell 25, 69-80. doi: 
10.1016/j.devcel.2013.03.002 

Wu, P., Zuo, X., Deng, EL, Liu, X., Liu, L., and Ji, A. (2013). Roles of long noncoding 
RNAs in brain development, functional diversification and neurodegenerative 
diseases. Brain Res. Bull 97, 69-80. doi: 10.1016/j.brainresbull.2013.06.001 

Wu, Q., Wang, Y, Cao, M., Pantaleo, V., Burgyan, J., Li, W. X., et al. (2012). 
Homology-independent discovery of replicating pathogenic circular RNAs by 
deep sequencing and a new computational algorithm. Proc. Natl. Acad. Sci. 
U.S.A. 109, 3938-3943. doi: 10.1073/pnas.lll7815109 

Yang, B., Lin, EL, Xiao, J., Lu, Y, Luo, X., Li, B., et al. (2007). The muscle-specific 
microRNA miR-1 regulates cardiac arrhythmogenic potential by targeting GJA1 
and KCNJ2. Nat. Med. 13, 486-491. doi: 10.1038/nml569 

Yekta, S., Shih, I. EL, and Bartel, D. P. (2004). Micro RNA-directed cleavage of 
HOXB8 mRNA. Science 304, 594-596. doi: 10.1126/science.l097434 

Zimonjic, D. B., Keck, C. L., Thorgeirsson, S. S., and Popescu, N. C. (1999). Novel 
recurrent genetic imbalances in human hepatocellular carcinoma cell lines 
identified by comparative genomic hybridization. Hepatology 29, 1208-1214. 
doi: 10.1002/hep.510290410 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 27 September 2013; accepted: 07 January 2014; published online: 30 January 
2014. 

Citation: Kartha RV and Subramanian S (2014) Competing endogenous RNAs 
(ceRNAs): new entrants to the intricacies of gene regulation. Front. Genet. 5:8. doi: 
10.3389/fgene.2014.00008 

This article was submitted to Non-Coding RNA, a section of the journal Frontiers in 
Genetics. 

Copyright © 2014 Kartha and Subramanian. This is an open-access article dis- 
tributed under the terms of the Creative Commons Attribution License (CC BY). The 
use, distribution or reproduction in other forums is permitted, provided the original 
author(s) or licensor are credited and that the original publication in this jour- 
nal is cited, in accordance with accepted academic practice. No use, distribution or 
reproduction is permitted which does not comply with these terms. 



www.frontiersin.org 



January 2014 | Volume 5 | Article 8 | 9 



